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Abstract 

ZEUS inclusive diffractive cross-section measurements have been used in a 
DGLAP next-to-leading-order QCD analysis to extract the diffractive parton 
distribution functions. Data on diffractive dijet production in deep inelastic 
scattering have also been included to constrain the gluon density. Predictions 
based on the extracted parton densities are compared to diffractive charm and 
dijet photoproduction data. 
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1 Introduction 



Many aspects of diffractive interactions can be described in the framework of quantum 
chromodynamics (QCD) as long as a hard scale is present, so that perturbative techniques 
can be used and the dynamics of the processes can be formulated in terms of quarks 
and gluons. HERA data have contributed significantly to the understanding of such in- 
teractions, since events characterised by the diffractive dissociation of virtual photons, 
7*p — > Xp, constitute a large fraction (« 10%) of the visible cross section in deep inelas- 
tic scattering (DIS). Diffractive reactions in DIS are a tool to investigate low- momentum 
partons in the proton, notably through the study of diffractive parton distribution func- 
tions (DPDFs). The latter are the densities of partons in the proton when the final state 
of the process contains a fast proton of specified four-momentum. A precise knowledge of 
the DPDFs is an essential input to predictions of hard diffractive processes at the LHC. 

Several recent sets of DPDFs [1-4] have been determined in global fits using the con- 
ventional DGLAP formalism [5-7] in next-to-leading-order (NLO) QCD. The diffractive 
structure function, F 2 D , which dominates the cross section, is directly sensitive to the 
quark density, whereas the gluon density is only indirectly constrained via scaling viola- 
tions in the inclusive diffractive DIS cross sections. The inclusion of diffractive dijet data 
provides an additional constraint on the gluon density, since gluons directly contribute to 
jet production through the boson-gluon fusion process [8]. 

For this paper, recently published ZEUS inclusive diffractive data [9] and diffractive DIS 
dijet data [10] were used to extract the DPDFs. The inclusive data were first fitted alone 
and then in combination with the dijet data. The results are compared to HI fits [1,2], as 
well as to ZEUS diffractive charm data [11] and to ZEUS diffractive dijet photoproduction 
data [12]. 

2 Theoretical framework 

The cross section for diffractive DIS, ep — > eXp, in the one-photon-exchange approxima- 
ted in terms of the di 

da ep ^ eXp 2na 2 



tion, can be expressed in terms of the diffractive reduced cross-section af^\ 



dpdQ 2 dx P (3Q A 



1 + (1-J/) 



<T?W(P,Cf,X P ), (1) 



which depends on the diffractive structure functions, as 



a^\(3,Q\xjp) = F^ 3 \{3, Q\ z P )- f iff (3) (/3, Q\ x F ). (2) 

i + (i - yy 



i 



The kinematic variables used in Eqs. ([ID and ([2D are defined as follows: 

• Q 2 = —q 2 = —(k — k') 2 , the negative invariant-mass squared of the exchanged virtual 
photon, where q = k — k' is the difference of the four-momenta of the incoming and 
outgoing leptons; 

• xp = (P — P') ■ q/P ■ q, the fraction of the momentum of the proton carried by 
the diffractive exchange, where P and P' are the four-momenta of the incoming and 
outgoing protons, respectively; 

• (3 = Q 2 /2(P — P') ■ q, the Bjorken variable defined for the diffractive exchange; 

• the inelasticity y — (q ■ P)/(k • P). 

The four-momentum exchanged at the proton vertex, \t\, is integrated over in Eqs. (JTJ) 
and (E]). 

The QCD factorisation theorem [13-16] allows the diffractive structure functions, F^f\ 
to be expressed in terms of convolutions of coefficient functions and DPDFs: 

1 

F$\(3,Q 2 ,xp) = J f C V^) fi D (*>**&)> ( 3 ) 
i P 

where the sum runs over partons of type i and z is the longitudinal momentum fraction of 
the parton entering the hard subprocess with respect to the diffractive exchange. In the 
lowest-order quark-parton model process, z = j3. The inclusion of higher-order processes 
leads to /3 < z. The coefficient functions CW,, are the same as in inclusive DIS. In 
analogy to the usual parton distribution functions, the DPDFs ff(z, xjp; Q 2 ) are densities 
of partons of type i with fractional momentum zxjp in a proton, probed with resolution Q 2 
in a process with a fast proton in the final state with fractional momentum (1 — xp). The 
dependence of the DPDFs on the scale Q 2 is given by the DGLAP evolution equations. 

Proton-vertex factorisation [17] was adopted to model the xp dependence of the DPDFs. 
Two contributions were assumed, called Pomeron and Reggeon, separately factorisable 
into a term depending only on xp and a term depending only on z and Q 2 , 

f t D (z, xp; Q 2 ) = f F (x P )fi(z, Q 2 ) + fnMf^iz, Q 2 ). (4) 

The flux-factors fp and fp describe the emission of the Pomeron and Reggeon from the 
proton. Such an assumption was shown [9] to work to a good approximation for the data 
used in this analysis. 
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3 Analysis method 



The DGLAP evolution equations yield the distributions fi(z, Q 2 ) of the quarks and gluons 
at all values of Q 2 , provided the DPDFs are parameterised as functions of z at some 
starting scale Qq. The input parameters were fitted to the data by minimising a x 2 
function [8]. Correlated systematic uncertainties were taken into account by using the 
method described in an earlier ZEUS publication [18]. 

The QCD evolution was performed with the programs Qcdnum [19] as well as a newly 
developed package, QcdcQ. The strong coupling constant was set to a s (M z ) = 0.118. 
The contribution from heavy quarks was treated within the general-mass variable-flavour- 
number scheme of Thorne and Roberts (TR-VFNS) [20], which interpolates between the 
threshold and the high-Q 2 behaviour for heavy quarks. The values of the heavy-quark 
masses used were m c = 1.35 GeV and mj = 4.3 GeV [8]. The influence of was 
accounted for through its NLO dependence on the parton densities. 

The NLO QCD predictions for the diffractive structure functions were obtained by con- 
voluting the DPDFs with the QCD coefficient functions. Predictions for the inclusive 
diffractive reduced cross sections were obtained using Q as the factorisation and renor- 
malisation scales; for the dijet cross sections, Q was taken as the factorisation scale and 
the transverse energy of the leading jet, -E 1 ^*, as the renormalisation scale. Predictions 
for the dijet cross sections at the parton level were performed at NLO with Disent [21] 
and NlojetH — \- [22]. The two programs agree within 5%. These programs can deal with 
an arbitrary number of flavours but treat all quarks as massless. Thus they match the 
TR-VFNS at scales much larger than the quark masses, but may give imprecise results 
for scales close to the mass thresholds. The predictions were corrected for hadronisation 
effects [23]. 

3.1 Parameterisation of the DPDFs 

The DPDFs were modelled at the starting scale Qq = 1.8 GeV 2 in terms of quark singlet, 
/+ = J2q(fq + fq)i an d gluon, f g , distributions. The neutrality of the diffractive exchange 
requires fg = f q for all flavours. The light-quark distributions were assumed to be equal, 
fu — fd — fs- The distributions of the c and b quarks were generated dynamically at the 
scale Q above the corresponding mass threshold, i.e. no intrinsic charm or bottom were 
assumed. At the starting scale, chosen to be below the charm threshold, the quark-singlet 
parameterisation is summed over the light-quark distributions, /+ = 6f q , where q denotes 



Computer code developed by W. Slominski. 
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any of u, d, s and their antiquarks. The distributions were parameterised at Ql as 

zfd >u ,s{z,Ql) = A q z B "{\- zf\ 

zf 9 (z,Ql)=A g z B *(l-zf*. (5) 

An additional factor, e~ , was included to ensure that the distributions vanish for 
z — > 1 even for negative values of C q>g . 

The xp dependence was parameterised using Pomeron and Reggeon fluxes 

Xjp 

with linear trajectories atp^t) = of ( .r(0) + a'jp K t, where t is the four-momentum transfer 
at the proton vertex. The flux factors in Eq. (J6j) were integrated over t between —1 GeV 2 
and the kinematically allowed maximum value, as for a previous ZEUS publication [9]. 
The values of the parameters which were fixed in the fits are summarised in Table [TJ 
The value of the normalisation parameter A$> was absorbed in A q , g . The Reggeon parton 
densities, ff 1 in Eq. ([41), were taken from a parameterisation derived from fits to pion 
structure-function data [24]. 

In total, 9 parameters were left free in the fits: A q>g , B qg , C qt9 , the Pomeron and Reggeon 
intercepts, ajp(O) and aj?(0), and the normalisation of the Reggeon term, Am- 



3.2 Data 

Inclusive diffractive data were fitted alone as well as in combination with a sample of 
diffractive dijet data. The inclusive diffractive reduced-cross-section data [9] used in the 
fits were selected using two methods: the requirement of a large rapidity gap between the 
final-state proton and the rest of the hadronic system (LRG sample) and the detection 
of the final-state proton in the ZEUS leading proton spectrometer (LPS sample j§. These 
data cover photon-proton centre-of-mass energies in the range 40 < W < 240 GeV, photon 
virtualities in the range 2 < Q 2 < 305 GeV 2 (LRG) and 2 < Q 2 < 120 GeV 2 (LPS) and 
hadronic-final-state masses in the range 2 < Mx < 25 GeV (LRG) and 2 < Mx < 40 GeV 
(LPS). They span the following X]P ranges: 0.0002 < x P < 0.02 (LRG) and 0.002 < 
xp < 0.1 (LPS). Both samples are corrected for background from double-diffractive events 
ep — > eXN, in which the proton also dissociates into a low-mass state, N. The extracted 

2 The ZEUS measurements obtained with a third selection method based on the analysis of the distri- 
bution of the hadronic-final-state mass (Mx method) [25,26] were not included in the fits as this data 
sample is highly statistically correlated with the LRG sample. 
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DPDFs hence correspond to the single-diffractive reaction ep — > eXp. Only inclusive 
data with Q 2 > 5 GeV 2 were included in the fits. In order to minimise overlap between 
the two samples, LPS data with xjp < 0.02 were excluded from the fit. Reduced cross 
sections in bins of Q 2 , xjp and (3 were used. The total number of fitted data points was 
229 from the LRG sample and 36 from the LPS sample. 

The diffractive dijet sample in DIS [10] used in the fits covers transverse energies Ej, > 
4 GeV and xp < 0.03. Differential cross sections as a function of z^ s , an estimator of z, 
in different Q 2 bins were used. The variable z°^ s was calculated , Tji , where 

Mjj is the invariant mass of the dijet system. This sample provided 28 additional points. 

The overlaps between the LRG and dijet samples and between the LPS and dijet samples 
are small. Therefore, the three data sets were considered statistically independent. The 
following sources of correlated systematic uncertainties were considered: 

• LRG data: energy scale, reweighting of the simulation in £p, variation of the en- 
ergy threshold on the most forward energy-flow object and proton-dissociation back- 
ground [9]; 

• LPS data: energy scale, reweighting of the simulation in xjp and t and proton- 
dissociation background [9]; 

• dijet data: energy scale and proton-dissociation background [10]. 

The energy scale was taken to be common for all the data sets. The normalisation uncer- 
tainty due to the proton-dissociation background was taken as fully correlated between 
the LRG and dijet samples. 



4 Results 

Fits were performed to the data sets described in Section 13.21 Two different parame- 
terisations of the gluon density at the starting scale were used. This is summarised in 
Table [2j 

Two fits were performed to the LRG+LPS data: 'Standard', with A g , B g and C g as free 
parameters; 'Constant', with B g = C g = 0. The corresponding sets of DPDFs are referred 
to as "ZEUS DPDF S" and "ZEUS DPDF C" , respectively. 

Only data with Q 2 > Q^ in = 5 GeV 2 could be fitted within the combined assumptions of 
DGLAP evolution and proton- vertex factorisation. The quality of the fit drops rapidly 
for Q^in < 5 GeV 2 , as shown in Fig. Q} where xV n df obtained for the fit ZEUS DPDF S, 
for statistical reasons restricted to the LRG data, is presented as a function of Q 2 lin . 
Although the proton-vertex factorisation assumption worked reasonably well at low Q 2 



5 



for the LRG data [9], the combination of this assumption with the DGLAP evolution in 
Q 2 breaks down below Q 2 of 5 GeV 2 . In contrast, in fully inclusive DIS, DGLAP fits were 
performed successfully [18] down to Q 2 of 2.5 GeV 2 . 

The results of the fits are listed in Table [31 In the minimisation procedure, the x 2 was 
calculated using only the statistical and uncorrelated systematic errors. The correlated 
systematic uncertainties were treated according to the "offset method" [18] and were in- 
cluded in the determination of the full experimental uncertainty. To evaluate the goodness 
of the fit, the \ 2 shown in Table [3] was recalculated by adding the statistical, uncorrelated 
and correlated systematic errors in quadrature. The number of degrees of freedom in this 
case corresponds to the number of fitted points. 

The inclusive data are sensitive to all three quark- related parameters, A q , B q and C q . 
However, they show very little sensitivity to the gluon shape (B g and C g ). This follows 
from being directly sensitive only to quarks. The values obtained for ajp(O), q;jr(0) 

and A]r in all fits are fully consistent with those extracted from a fit [9] to the same data 
based on the proton- vertex factorisation assumption of Eq. (jlj). 

Fits ZEUS DPDF C and S are of equally good quality and the predicted reduced cross 
sections are indistinguishable. Figures [2HH show fit S compared to the LRG and LPS data. 
Both data samples are well described. For Q 2 < 5 GeV 2 , the predictions are extrapolated 
and underestimate the LRG as well as the LPS data for xp < 0.005. The fit is above the 
LPS data in the \ow-(3 region, where there are no LRG data. 

The quark and gluon densities, zf q (z, Q 2 ) and zf g (z, Q 2 ), from fit S are shown with their 
experimental uncertainties in the upper (lower) part of Fig. [5] for Q 2 = 6, 20, 60 and 200 
GeV 2 . The DPDFs from fit C are also shown. The relative normalisation of the Pomeron 
flux, fjp, and of the DPDFs, fi(z, Q 2 ), in Eq. (j4j) is arbitrary; for Fig. [5]the normalisation 
was chosen such that at xjp = 0.003 the quantity £p/jp(xjp) = 1. 

The quark distributions are very similar for the two fits while the gluon densities are 
significantly different. Gluons from fit S grow rapidly at high z, while those from fit C 
vanish as z — > 1 in a smoother way. The large discrepancy demonstrates the low sensitivity 
of the inclusive data to gluons. To constrain the gluons better, a more exclusive process 
is needed where photon-gluon fusion contributes at leading order. 

Predictions based on ZEUS DPDF S and C are compared to the double-differential diffrac- 
tive dijet cross section [10] in Fig. [61 The uncertainty due to the variation of the renor- 
malisation scale in the calculation between O.BE 3 ^ and 2E^ t is also shown. At high z°p s 
the predictions based on fit C give a good description of the dijet data throughout the 
whole kinematic region. This is not the case for fit S. Figure [6] clearly shows that the dijet 
data are sensitive to the gluon density. 

A third fit was performed to the LRG+LPS data in combination with the dijet data. The 
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corresponding DPDFs are referred to as "ZEUS DPDF SJ" and the resulting parameters 
are also given in Table [3j This fit is indistinguishable from fit S when compared to the 
inclusive data and also provides a remarkably good description of the dijet data, as shown 
in Fig. 

The quark and gluon densities from fit SJ are shown with their experimental uncertainties 
in Fig. [8]for Q 2 = 6, 20, 60 and 200 GeV 2 . The result of fit C is also shown for comparison. 
Again, the plotted quantities are zf q:9 (z,Q 2 ) with the normalisation xjpfjp(x]p) = 1 at 
xjp = 0.003. The decrease in the uncertainty on the gluon distribution with respect to 
the fit without jet data (Fig. [5]) is clearly seen. Combining the inclusive and dijet data 
constrains the gluon and the quark densities with a comparable precision across the whole 
z range. 

Figure [9] shows the Q 2 dependence of the fraction of the longitudinal momentum of the 
diffractive exchange carried by the gluons, g{ rajC , according to fit ZEUS DPDF SJ, inte- 
grated over the range 10~ 5 < z < 1. Over the wide range 5 < Q 2 < 300 GeV 2 , the 
fraction amounts to approximately 60%. The fall with Q 2 is a direct consequence of the 
DGLAP evolution which forces gf rac to approach « 0.5 at high Q 2 , while the slope change 
at Q = mj reflects the change in the number of flavours from four to five. 

4.1 DPDF uncertainties 

The following sources of uncertainties were investigated: 

• the starting scale Qq. The value 1.8 GeV 2 was chosen as it minimises the x 2 - Varying 
Ql for fit C between 1.6 and 2 GeV 2 yielded a \ 2 between 1.18 and 1.20; the DPDFs 
did not change significantly; 

• the fixed parameters in the fits (Tabled]). Variations within the measurement errors 
resulted in a simple scaling of the fluxes integrated over t, absorbed into the normali- 
sation parameters A q , A g and Ajr, with negligible effect on the DPDFs; 

• the renormalisation scale dependence. The scale /Ir for dijet data was taken as O.dE 3 ^ 
and 2E 3 r£ t , whereas it was kept as Q for the inclusive data. The effect on the parton 
densities was within 5% for light quarks, 15% for c and b and 30% for gluons, while 
the x 2 increased significantly; 

• the masses of the charm and beauty quarks. The nominal values of m c = 1.35 GeV and 
rrib = 4.3 GeV were varied in the ranges 1.35 < m c < 1.75 GeV and 4.3 < m& < 5 GeV. 
Neither the quark nor gluon distributions were sensitive to variations of wij, whereas 
m c produced an effect comparable to the experimental uncertainty. The \ 2 value 
changed only slightly, reaching the minimum at the nominal mass values. 
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These theoretical uncertainties are not shown in Figs. [S] and O Given are the total 
experimental uncertainties as determined with the offset method [18]. 

4.2 Comparison to other HERA DPDFs 

An earlier ZEUS DPDF fit to inclusive diffractive measurements combined with data 
on charm production [27] did not include data points with values of (3 > 0.4. Also, in 
the current analysis a far larger data sample is used. This previous analysis is therefore 
superseded by the results of the present paper. 

Several DPDF fits are available from the HI Collaboration, among which the most recent 
are based on an inclusive diffractive sample [1] and on a combination of data from inclusive 
and diffractive dijets in DIS [2]. In Fig. [10l the diffractive reduced cross sections from 
the ZEUS LRG sample as a function of Q 2 for fixed (3 and xp are compared to the fit 
ZEUS DPDF SJ as well as to the predictions from "HI 2006 Fit B" [1]. The latter, 
extracted for masses of the dissociative system < 1.6 GeV, were multiplied by the 
scaling factor 0.81 [1] to account for the correction to the ZEUS elastic case Mn = m p . 
For Q 2 < 5 GeV 2 in the ZEUS case and for Q 2 < 8.5 GeV 2 in the HI case, the fits are 
extrapolated. For /3 < 0.2, the two fits agree in shape throughout the fitted range, but 
the normalisation of the ZEUS curve is above that of HI. At higher j3 and where the 
predictions are extrapolated the agreement worsens. These features reflect the degree of 
consistency between the ZEUS and HI data [28]. 

4.3 Comparison to charm and dijet photoproduction data 

Predictions from the fit ZEUS DPDF SJ are compared in Fig. [11] to data [11] on the 
charm contribution to the diffractive structure function multiplied by xjp, X]pF®^ cc , for 
xjp of 0.004 and 0.02 and Q 2 of 4 and 25 GeV 2 . The predictions are in fair agreement 
with the data. 

Figures fT2h and fT3h show the prediction from the fit ZEUS DPDF S J compared to the 
diffractive dijet photoproduction data [12] as a function of x° bs , the fraction of the pho- 
ton energy invested in producing the dijet system, and of the tranverse energy Eif? of 
the leading jet, respectively. Dijet photoproduction at leading order in QCD proceeds 
through two type of processes: direct processes (x° bs ~ 1), in which the exchanged pho- 
ton interacts as a point-like particle with the partons from the diffractive exchange; and 
resolved processes (x° bs < 1), in which the photon behaves as a source of partons which 
interact with the diffractive exchange. For the latter, as in hadron-hadron interactions, 
QCD factorisation is not expected to hold [13]. Comparing calculations based on DPDFs 
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extracted from other diffractive processes to diffractive dijet photoproduction data pro- 
vides a valuable test of QCD factorisation. The predictions shown in Figs. fT2a and [T3h . 
obtained with the program of Klasen and Kramer [29], agree adequately with the data 
over the whole x° bs and Ej? ranged. A further confirmation is provided by Figs. [T2b and 
[13b, where the ratio of the data and the predictions is presented as a function of x° bs 
and Ej?, respectively. The ratio is consistent with unity. This reinforces the conclusion 
of an earlier ZEUS publication [12], where the data were found to be compatible with no 
suppression either of the resolved component, or of both components globally. The slight 
normalisation difference between the predictions shown here and those of the earlier ZEUS 
publication [12] is due to the usage of different DPDFs ("HI 2006 Fit B"), extracted in a 
fixed-flavour-number scheme. 

Also shown in Figs. fT2l and [TBI are the predictions from "HI Fit 2007 Jets" [2], multiplied 
by the scaling factor 0.81 [1] to account for the correction to = m p . They overlap 
with the ZEUS predictions in the lowest and highest x° bs bins. At high E^ and low 
x° bs , the two predictions differ by as much as 20%. 

5 Summary and conclusions 

ZEUS diffractive inclusive cross-section data, together with data on diffractive dijet pro- 
duction in DIS, have been used in a NLO DGLAP QCD analysis to determine the diffrac- 
tive parton distribution functions. Only data with Q 2 > 5 GeV 2 could be fitted within the 
combined framework of DGLAP evolution and proton- vertex factorisation. The extracted 
DPDFs correspond to the single-diffractive reaction with a proton in the final state and 
are valid in the region \t\ < 1 GeV 2 , M x > 2 GeV, x F < 0.1. 

NLO QCD predictions based on diffractive parton densities extracted from the inclusive 
data provide a good, simultaneous description of the inclusive and dijet data. The quark 
densities are well constrained by the inclusive data, whereas the sensitivity to the gluon 
density decreases as the momentum fraction z increases. 

The inclusion of the dijet data provides an additional constraint on the gluons, allowing 
the determination of both the quark and gluon densities with good accuracy. 

Predictions based on the extracted parton densities give a fair description of diffractive 
charm production data; also, they are in good agreement with diffractive dijet photopro- 
duction cross sections over the whole kinematic region, thus indicating no suppression of 
the resolved component. 



3 The program of Frixione and Ridolfi [30] was also used and gave the same results. 
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Parameter 


Fixed to 


Measurement 


Ref. 


a'jp 

Bp 
Bm 




0.9 GeV" 2 
7.0 GeV~ 2 
2.0 GeV" 2 


-0.01 ± 0.06(stat.)+°;^(syst.) ± 0.04(model) GeV" 2 
0.90 ± 0.10 GeV" 2 
7.1 ± 0.7(stat.)^(syst.) GeV^ 2 
2.0±2.0GeV~ 2 


[9] 
[31] 

[9] 
[31] 



Table 1: The values of the parameters fixed in the fits and the measurements 
providing this input. 



Fit name 


Data set 


zg{z) 


ZEUS DPDF S 
ZEUS DPDF C 
ZEUS DPDF SJ 


LRG + LPS 
LRG + LPS 

LRG + LPS + DIS dijets 


A g z Ba (l - z)° a 

Ag Z B °(l ~ Z) C9 



Table 2: Data sets used and parameterisation of the gluon density at the starting 
scale for the different fits. 



Parameter 


Fit value 


Fit value 


Fit value 




DPDF S 




DPDF C 


DPDF SJ 


A i 


0.135 


± 


0.025 


0.161 


± 


0.030 


0.151 


± 


0.020 


B q 


1.34 


± 


0.05 


1.25 


± 


0.03 


1.23 


± 


0.04 


c q 


0.340 


± 


0.043 


0.358 


± 


0.043 


0.332 


± 


0.049 




0.131 


± 


0.035 


0.434 


± 


0.074 


0.301 


± 


0.025 




-0.422 


± 


0.066 









-0.161 


± 


0.051 


c 9 


-0.725 


± 


0.082 









-0.232 


± 


0.058 


«7P(0) 


1.12 


± 


0.02 


1.11 


± 


0.02 


1.11 


± 


0.02 


«k(0) 


0.732 


± 


0.031 


0.668 


± 


0.040 


0.699 


± 


0.043 


A m 


2.50 


± 


0.52 


3.41 


± 


1.27 


2.70 


± 


0.66 


X 2 /ndf 


315/265 




1.19 


312/265 




1.18 


336/293 




1.15 



Table 3: Parameters obtained with the different fits and their experimental un- 
certainties. 
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Figure 1: Dependence of x 1 /ndf for the fit ZEUS DPDF S on the value of the 
minimum Q 2 , Qmm> of the LRG data entering the fit. 
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Figure 2: The fit ZEUS DPDF S compared to the ZEUS LRG data [9] as a 
function of xp for different (3 and Q 2 values at low Q 2 . Where visible, the inner 
error bars show the statistical uncertainties and the full bars indicate the statistical 
and systematic uncertainties added in quadrature. The dashed lines represent the 
DGLAP extrapolation beyond the fitted region. 
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Figure 3: The fit ZEUS DPDF S compared to the ZEUS LRG data [9] as a 
function of xp for different (3 and Q 2 values at high Q 2 . Where visible, the inner 
error bars show the statistical uncertainties and the full bars indicate the statistical 
and systematic uncertainties added in quadrature. 
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Figure 4: The fit ZEUS DPDF S compared to the ZEUS LPS data [9] as a function 
of Xp for different (3 and Q 2 values. The inner error bars show the statistical 
uncertainties and the full bars indicate the statistical and systematic uncertainties 
added in quadrature. The dashed lines represent the DGLAP extrapolation beyond 
the fitted region. 
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Figure 5: Four upper (lower) plots: the quark (gluon) distributions obtained from 
fits ZEUS DPDF S (continuous line) and ZEUS DPDF C (dashed line), shown 
for four different values of Q 2 . The shaded error bands show the experimental 
uncertainty. 
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Figure 6: Predictions based on ZEUS DPDF S ( continuous line) and ZEUS DPDF 
C (dashed line) compared to the ZEUS diffractive dijet data [10] as a function of 
z ip s f or different Q 2 values. The inner error bars show the statistical uncertain- 
ties and the full bars indicate the statistical and systematic uncertainties added 
in quadrature. The dark shaded bands indicate the jet energy scale uncertainty. 
The light shaded bands at the bottom of each plot show the renormalisation scale 
uncertainty. 
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Figure 7: JTie fit ZEUS DPDF SJ compared to the ZEUS diffractive dijet data [10] 
as a function of z°^ s for different Q 2 values. The inner error bars show the sta- 
tistical uncertainties and the full bars indicate the statistical and systematic uncer- 
tainties added in quadrature. The dark shaded bands indicate the jet energy scale 
uncertainty. The light shaded bands show the DPDF experimental uncertainties. 
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Figure 8: Four upper (lower) plots: the quark (gluon) distributions obtained from 
fits ZEUS DPDF SJ (continuous line) and ZEUS DPDF C (dashed line), shown 
for four different values of Q 2 . The shaded error bands show the experimental 
uncertainty. 
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Figure 10: The fit ZEUS DPDF SJ (continuous line) and HI Fit B [1] (dashed- 
dotted line) compared to the ZEUS LRG data [9] at (a) xp = 0.003 and (b) xp = 
0.01 as a function of Q 2 for different j3 values. Where visible, the inner error 
bars show the statistical uncertainties and the full bars indicate the statistical and 
systematic uncertainties added in quadrature. The HI predictions are corrected 
to Mn = m p via the scaling factor 0.81. The dashed (dotted) lines represent the 
DGLAP extrapolation beyond the ZEUS (HI) fitted region. 
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Figure 11: ZEUS DPDF SJ predictions compared to the ZEUS measurement 
of the charm contribution to the diffractive structure function multiplied by xjp, 
xjpF2^ cc [11] as a function of (3 for different Q 2 and xjp values. The inner error 
bars show the statistical uncertainties and the full bars indicate the statistical and 
systematic uncertainties added in quadrature. The dashed lines represent the ex- 
trapolation of the predictions beyond the kinematic region where they were obtained. 
The shaded bands show the DPDF experimental uncertainty. 
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Figure 12: (a) ZEUS DPDF S J predictions compared to the ZEUS diffractive 
dijet photoproduction data [12] as a function of a;° bs . The predictions from HI Fit 
2007 Jets [2] are also shown, corrected to = m p via the scaling factor 0.81. 
(b) Ratio of the data and of the HI predictions to the ZEUS DPDF SJ predictions. 
The inner error bars show the statistical uncertainties and the full bars indicate 
the statistical and systematic uncertainties added in quadrature. The dark shaded 
bands indicate the jet energy scale uncertainty. The light shaded bands show the 
DPDF experimental uncertainty. The shaded band at the bottom of (a) shows the 
renormalisation scale uncertainty. 
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Figure 13: (a) ZEUS DPDF S J predictions compared to the ZEUS diffractive 
dijet photoproduction data [12] as a function of ■ The predictions from HI Fit 
2007 Jets [2] are also shown, corrected to = m p via the scaling factor 0.81. (b) 
Ratio of the data and of the HI predictions to ZEUS DPDF SJ predictions. Other 
details as in caption to Fig. [7H 
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